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FOREWORD 


The  work  documented  in  this  report  is  part  of  a  study  of 
the  travel-time  and  amplitude  anomalies  observed  at  the  Large 
Aperture  Seismic  Array  (LASA)  in  Montana.  This  report  super¬ 
cedes  Report  No.  LL-4  in  that  additional  data  and  data  quality 
control  are  used  in  this  report. 

The  work  was  performed  by  the  Applied  Research  Section, 

Earth  Sciences,  a  Teledyne  Company,  314  Montgomery  Street, 
Alexandria,  Virginia,  under  Lincoln  Laboratory  Contract  No.BB-246 
This  report  was  written  by  P.  W.  Broome,  F.  A-  Klappenberger 
Do  Eo  Frankowskio  Assistance  was  provided  by  A.  L.  Kurtz,  V.  R. 
McLamore,  R.  C.  Mierley,  P.  Ao  Santiago,  H.  N.  Johnson,  and  F.  J. 
Whited  o  Dr.  E.  F.  Chibuns  served  as  a  consultant.  The  project 
director  was  P.  W.  Broome. 


ABSTRACT 


Average  amplitude  anomalies  are  presented  as  functions  of 
source  region.  Data  from  approximately  300  events  were  used  in 
this  analysis. 


Accepted  for  the  Air  Force 

Franklin  C.  Hudson 

Chief,  Lincoln  Laboratory  Office 


INTRODUCTION 


1  . 

Observed  short  period  teleseismic  P-phase  amplitudes 
vary  greatly,  even  between  seismometers  that  are  relatively 
closely  spaced.  These  variations  exceed  anything  that  can  be 
easily  attributed  to  instrumentation,  radiation  pattern  ef¬ 
fects,  or  obvious  geologic  factors.  These  variations  have 
been  termed  "amplitude  anomalies"  By  observing  a  large 
number  of  events,  certain  consistencies  are  found,  in  par¬ 
ticular,  some  stations  exhibit  a  preference  for  events  from 
a  given  source  region 

The  object  of  this  report  is  to  define  the  LASA  ampli¬ 
tude  anomalies  as  functions  of  distance  and  azimuth  to  the 
source  region.  Statistical  qualifications  on  the  observed 
anomalies  and  on  the  data  from  approximately  300  events 
(Figure  1)  which  were  used  to  derive  these  anomalies  are  also 
given  . 

2  PROCEDURE 

2 . 1  Data  Collection 

Figure  2  shows  a  typical  record  of  an  event  as 
seen  on  the  LASA  16mm  Develocorder  film.  Two  such  films  re¬ 
cord  the  earth  motion  at  the  center  seismometer  of  each  of 
the  21  LASA  subarrays.  The  maximum  P-wave  amplitude  and  period 
within  the  first  3  or  4  seconds  of  the  signal  were  measured 
and  converted  to  ground  motion* .  Since  all  sensors  do  not  re- 

-  1  - 


*The  maximum  P-wave  amplitude  is  defined  to  be  that  phase  which 
has  the  largest  amplitude,  on  the  average,  for  all  recording 
seismometers 


Figure  1.  Event  Location  Chart  with  Distance-Azimuth  Grid 
(Distance  in  Kilometers,  Azimuth  in  Degrees) 
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F igure  2 . 


Typical  Record  Showing  Amplitude  Measurements 


cord  with  the  same  signal-to-noise  characteristics,  each  trace 
was  given  a  quality  grade  by  the  analyst  when  read. 

Quality  Grade  1  =  a  perfect  trace. 

Quality  Grade  2  =  a  trace  on  which  the  signal  is  easily 

read  and  has  a  similar  shape  (signature) 
to  the  other  traces. 

Quality  Grade  3  =  a  trace  on  which  the  signal  is  seen 

but  it  is  not  distinct  and/or  has  a 
markedly  dissimilar  signature  to  the 
other  traces. 

Quality  Grade  4  =  a  trace  on  which  the  signal  can  not  be 

seen . 

In  this  report  all  but  Grade  4  data  were  used. 

These  data  were  then  keypunched  on  IBM  cards  as  shown  in 
Figure  3  for  use  in  further  processings. 

2.2  Data  Reduction 


In  order  to  compare  data  from  different  events, 

it  was  first  necessary  to  remove  the  effect  of  the  sources  on 

the  received  amplitudes.  In  doing  this,  we  have  assumed  that 

the  aperture  of  LASA  is  small  compared  to  any  lobe  structure 

in  the  radiation  pattern,  and  thus  that  all  seismometers  were 

sampling  the  same  radiation  field.  The  normalization  method 

which  we  used  was  to  normalize  the  amplitude  at  each  station 

for  each  event  to  the  geometric  mean  of  the  recorded  amplitudes 

for  that  event  Thus  if  Y  and  L. .  are  the  normalized  and 

13  13 

observed  amplitudes  for  event  i  and  station  j  respectively  and 
N  is  the  total  number  of  observed  amplitudes  for  event  i,  then 


Y 

13 


N 

r  n 

\j=1 


L  . 
13 


2 


a  an 


The  derivation  and  rationalization  of  this  normalization  method 
is  given  in  Appendix  A. 

A  computer  program  was  used  to  normalize  amplitudes 
for  each  event* .  A  sample  output  from  this  program  is  shown 
in  Figure  4.  (Additional  data  computed  in  the  program  are  also 
shown . ) 

2 . 3  Data  Display 

The  events  were  grouped  into  cells  by  distance  and 
azimuth.  With  the  assumption  that  the  amplitude  anomaly 
within  each  cell  is  constant,  a  computer  program  was  used  to 
compute  the  geometric  mean  of  the  normalized  amplitudes  for 
each  station  ** .  The  program  also  computed  a  measure  of  the 
deviation  about  the  estimated  mean  as  well  as  the  number  of 
events  recorded  for  each  station  and  cell***. 

Figure  5  shows  a  sample  of  this  information  as  the 
computer  output  listing  with  appropriate  descriptions  added. 

The  data  for  all  cells  and  subarrays  are  presented  in  Figures  6 
through  26 

3 .  ERROR  DETECTION 

Using  the  assumption  that  the  true  amplitude  anomalies 
are  constant  within  geographic  cells,  it  is  possible  to  ex¬ 
amine  the  scatter  of  the  normalized  amplitudes  in  each  cell 
for  the  purpose  of  detecting  errors  due  to  blunders  in  data 
handling  (and  other  possible  sources  of  error) .  To  detect 
and  eliminate  the  scatter  due  to  these  human  errors,  the  error 
model  study  described  in  Appendix  D  was  undertaken.  The  re¬ 
sults  of  this  study  indicate  that  the  signal  model  is  essential¬ 
ly  of  the  form  L  =  S  g.e. .,  where  L.  is  the  observed 

13  1^3  13  13 

-  3  - 

A  description  of  program  "ANOMALY"  is  given  in  Appendix  B 

**  Appendix  A 

***  A  description  of  program  "DISPLAY"  is  given  in  Appendix  C 
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Figure  4.  Sample  Output  from  Program  "ANOMALY 
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Figure  5.  Sample  Output  of  Program  "DISPLAY 
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Figure  15.  Subarray  D2 
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amplitude  for  event  1  and  station  j ,  S  is  the  radiated  signal 
of  event  i,  g  is  the  gain  factor  of  station  j  and  e  is  an 

3  13 

error  term  which  has  a  log-normal  distribution 

The  procedure  used  to  detect  errors  was  as  follows,  The 

2 

estimates  of  the  variances,  S  ,  computed  for  each  geographic 

cell  and  station  were  plotted  on  a  histogram  as  a  function  of 

the  number  of  events,  N,  in  that  cello  By  pooling  these  var- 

~2 

lances  a  very  good  estimate  S  ,  of  the  variance  of  the  error 

2  2 
population,  a  ,  was  found  Confidence  intervals  for  the  S  , 

2 

were  then  computed  and  all  S  which  fell  outside  of  the  upper 
95%  confidence  limit  were  questioned .  The  amplitude  measure¬ 
ments  involved  in  computing  these  outlier-estimates  were  re¬ 
examined  for  blunders.  (Readings  which  did  not  clearly  include 
blunders  were  not  changed,  so  this  re-examination  of  the  data 
should  not  bias  the  final  results  )  The  percentage  of  reading 
blunders  so  detected  was  less  than  1%  of  the  total  readings 
taken.  With  these  errors  eliminated,  the  entire  process  was 
repeated  with  the  results  shown  by  Figure  27.  On  this  second 
pass  the  total  number  of  outliers  above  the  upper  confidence 
limit  is  1.1%  of  the  total  number  of  estimates. 


4.  CONCLUSIONS 

It  is  clear  from  the  results  of  this  study  that  certain 
stations  at  the  Large  Aperture  Seismic  Array  show  a  preference 
"on  the  average"  for  events  from  particular  regions.  Due  to 
the  high  variability  of  the  data,  however,  one  might  hypothe¬ 
size  that  the  averaging  yields  only  the  low  frequency  components 
of  an  anomaly  versus  distance-azimuth  function  which  contains 
in  addition,  a  considerable  amount  of  high  frequency  components. 
The  amount  of  available  data  is  not  adequate  to  provide  an 
answer  to  this  hypothesis.  It  is  interesting  to  note,  however. 


4 


that  the  average  anomalies  for  the  eight  Fiji  events,  used  in 
Appendix  A  to  establish  the  distribution  of  the  errors  have 
the  same  distribution  function  as  the  errors;  see  Figures  A3 
and  A4 . 

The  average  anomalies  are  large  in  many  cases,  with 
factors  of  3  to  1  between  stations  being  common.  As  such  they 
must  play  a  part  in  studies  involving  the  estimation  of  source 
radiation  patterns.  The  variability  of  the  data  might  also 
put  limits  on  the  degree  to  which  one  can  hope  to  estimate  these 
radiation  patterns  using  any  seismic  net:  a  trade-off  between 
resolution  and  accuracy  likely  exists. 

The  signal  model  which  was  derived  for  error  detection 
served  well  in  that  function,  however,  the  model  is  also  interest¬ 
ing  in  its  own  right  The  physical  significance  of  the  multipli¬ 
cative  error  factor  was  not  to  be  part  of  this  study,  and  for  this 
reason  no  effort  was  made  to  explain  it.  It  is  possible  that  the 
effect  could  have  some  bearing  upon  the  placement  of  new  arrays. 

The  data,  as  presented,  is  adequate  to  describe  the 
average  amplitude  anomalies  in  the  more  seismically  active 
regions  . 
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S  (standard  deviation) 


Figure  27. 


Distribution  of  Standard  Deviation  Estimates 
Computed  from  N  Event 


APPENDIX  A 


A .  1  INTRODUCTION 

The  signal  levels  recorded  at  the  stations,  aside  from 
statistical  variations,  will  be 


L  =  S  g 
il  il 


(A-l) 


where  is  the  signal  radiated  by  the  ith  event  and  g_.  is 
a  gain  associated  with  the  jth  station.  The  set  of  ratios 
of  these  levels  for  a  single  event  are  termed  "amplitude 
anomalies"  .  We  define  an  individual  anomaly  for  stations  -f 
and  m  as 


l mi 

note  that 
Am-M 


L(./L  . 
mi 


(Afmi } 


-1 


( A-2  ) 


(A-3) 


and  denote  the  set  of  these  as 


A  =  {a  .}  2_\  (A-4 ) 

l  l  -fmiJ  t  ,m-l 

All  elements  in  A^  are  equally  important  but  only  20 
of  the  441  are  independent.  Any  20  which  contain  each  station 
index  at  least  once  are  sufficient  to  construct  the  entire 
set  . 

In  analyzing  the  recorded  signals,  it  was  found  that 
A^  was  highly  variable.  To  reduce  this  variability,  averages 
were  taken  over  events  from  the  same  region.  This  pooling  of 


A-l 


events  is  possible  under  the  assumption  that  amplitude  anoma¬ 
lies  are  a  slowly  varying  function  of  distance  and  azimuth.* 
Thus,  if  the  region  is  sufficiently  small,  this  is  equivalent 
to  assuming  the  anomalies  are  constant  for  each  region. 

Several  approaches  were  taken  to  this  pooling  problem. 

Each  of  these  was  equivalent  in  the  sense  that  they  all  pro¬ 
duce  the  same  final  results.  A  summary  of  these  will  be 
presented  in  the  remainder  of  this  section. 

A  .  2 . 1  Method  1 

From  Equation  (A-3)  we  note  that  the  apparent 
size  of  an  anomaly  depends  upon  the  reference  station.  To 
eliminate  this  effect,  a  logarithmic  transformation  was  made. 

Since,  on  a  logarithmic  scale,  the  reference 
station  which  is  chosen  serves  only  to  translate  all  anoma¬ 
lies  equally,  we  can  dispense  with  a  reference  station  momen¬ 
tarily  and  concentrate  instead  on  a  set  of  translating  factors. 

First,  consider  the  problem  of  fitting  the  ob¬ 
served  amplitudes  for  a  single  event  to  a  selected  level  S. 

We  seek  a  multiplier  or  normalizing  factor,  G.  to  minimize 
the  total  deviation  from  S, 

v  2 

e.  -  )  (log  S-log  G,L  .)  (A-5) 

l  Li  l  lj 

j 

stat ions 


A-2 


■*"This  assumption  is  open  to  question.  It  is  clear  from  the 
results  of  this  study,  however,  that  certain  stations  do  show 
a  preference  ''on  the  average"  for  events  from  particular 
regions.  From  the  data,  though,  one  might  hypothesize  that 
the  averages  yield  only  the  slowly  varying  components  of  a 
distance-azimuth  function  which  contains  a  considerable  amount  of 
high  frequency  components.  The  amount  of  available  data  is  not 
adequate  to  answer  this  question. 


and  find  that 


G 

1 


=  S/ 


J 


TT 

3  =  1 


V 


a/j 


(A-6) 


where  J  is  the  total  number  of  stations. 

Second,  consider  minimizing  the  total  variation  between 
events.  We  define  this  total  variation  as 


L 


V' 

v  =  L 

j  l;k 

station  events 


(log  G  L 

1  13 


109  SLk  j  > 


(A-7  ) 


where  the  G  are  multiplicative  factors  as  before. 

This  minimum  will  be  found  among  the  solutions  to 
the  equations 


dV 

*G 


n 


0, 


n  1  ■  2  ,  o  •  • 


(A-8 ) 


or 


V  2 

L  G 

,  n 


(log  G  L  -log  G,  L.  .  ) 
n  n]  k  kj 


0 


(A-9) 


which  reduces  to 


I  log  G  J 
n 


Z  iog 

j 


n3 


(A— 10 ) 


where  I  and  J  are  the  total  number  of  events  and  stations 
respectively. 

By  setting 


Z  109  Gk  Lk: 

3  -k 


IJ  log  S , 


(A— 11 ) 


A- 3 


a  constant,  we  obtain 


J 

G  =  S/'  (  tt 

1  j-1 


L  ) 
ij 


V  1/J 


(A-12 ) 


as  before 


The  constant  S  can  be  chosen  arbitrarily.  The  weights 
G_^  are  independent  and  the  normalizing  process  is  equivalent  to 
fitting,  in  a  least-square  sense,  each  set  of  measurements,  L  , 
to  a  reference  level  Log  S  This  permits  addition  or  omission 
of  events  from  a  collection  without  disturbing  the  normalization 
factors  on  the  remaining  data.  In  this  way  we  obtain  a  set  of 
normalizing  factors  which  can  be  used  on  an  event  by  event  basis 
The  normalized  amplitudes  can  be  plotted  as  ordinates  on 
semi-log  graph  paper  with  stations  as  abscissa  on  the  linear 
axis.  The  resultant  plots  will  be  adjusted  in  a  least-squares 
sense  to  minimize  the  total  variation. 

To  use  these  data  for  computing  anomalies  we  first 
define  the  average  normalized  amplitude  at  each  station  as 

I 


log  g 


=  <  log  G  L  >  =  \  / 

D  i  iD  1  L 


log  G  L 

i  1 D 


(A-13) 


i-l 

events 


where  the  G  are  previously  determined  normalizing  factors.  Then 
the  average  anomaly  is  given  as 


< 


m£i 


log  G  L  > 
l  lm 


I 

TT 

i=l 


L  1/  1 
im 


<  log  G  L  . 

l  i/. 


> 


J 


(A-14) 


A -4 


which  is  the  geometric  mean  of  the  individual  anomaly  esti¬ 
mates.  This  reasonable  result,  along  with  the  desirable 
features  of  the  logarithmic  transformation  are  reason  enough 
for  using  this  method. 

A. 2. 2  Method  2 

Up  to  this  point  in  the  analysis  it  was  not  neces¬ 
sary  to  incorporate  an  error  model.  To  lend  further  justification 
to  the  above  method,  however,  it  was  felt  desirable  to  do  this. 

We  therefore  assumed  an  error  model  of  the  form 


L..  =  S  .  g.  e..  +  g  . 
ID  1  D  ID  ID 


(A-l 5 ) 


where  g. .  is  due  to  the  additive  factors  of  noise  and  reading 
ID 

errors  and  e .  accounts  for  variation  in  system  gain  due  to 
equipment  changes,  effects  of  wave-shapes,  etc. 

The  steps  taken  will  be  to  assume  and  then 
justify  two  hypotheses  which  were  put  forth  after  looking  at 
a  large  amount  of  data.  These  assumptions  are  that  the  addi¬ 
tive  noise  e  .  is  negligible  for  large  events  and  that  e_  has 
a  log  normal  distribution. 

By  neglecting  the  additive  error,  one  may 
linearize  the  model  to 


log  L  =  log  S  +  log  g .  +  log  e  . 
ID  i  D  ID 


(A-l 6 ) 


and  by  minimizing  the  sum  of  squares 


=  Y  (log  e.  ,)2  -  Y  (log  -g±l-  )2 


1/  D 


1 '  j 


S  .g . 
id 


(A-17) 


estimate  the  values  of  and  g_.  (albeit  not  independently)  . 


A- 5 


These  critical  values  of  gain  and  signal  level 
are  in  the  solutions  to  the  equations 


9E  _  ,  9E  _ 

-  0  and  — r —  -  0 


dS 

n 


9g 


m 


or 


V 


>  Log  (L  ,/S  g  .  )  =  0 
L>  nj  n  j 

1 


and 


V  Log  (L.  /S  g  )  =  0 
L  lm  l  m 


which  can  be  rewritten  as 


S  =  (  tt  L  ,/g  ) 


n 


1/J 


1=1 


nl  1 


and 


g  =  (  tt  L  /S .  ) 


m 


1/1 


i=l 


lm  i 


n  1 , 2  ,  .  .  .  ,  I , 

m  1,2,..., J 


( A-18 ) 


(A-19.1) 


(A-19.2) 


(A-20.1) 


(A-20.2) 


Then  by  adding  one  constraint,  namely  that 


J 

rr  g  .  =  S*" 


1=1 
we  have 


.  L  ,1/1 

(S  tt  im) 

i=l 
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IJ 

(  tt  L  .  ) 

ij-i  13 


1/IJ 


(A-21) 


(A-22) 


which  is  equivalent  to  Equation  (A-12)  when  1=1  and 

g  =  G  L 
m  l  un 


A-6 


The  anomalies  are,  as  before,  ratios  of  the 
estimates  of  the  gam  factors, 


<  A 


ml  j 


>  = 


g/g, 

m  L 


I 

tt 

i=l 


(L  /L  ) 
lm  i/. 


1/1 


(A-2  3 } 


(A-24 , 
A-14 ) 


A  3  Distribution  of  Errors 


To  establish  the  distribution  function  for  e. eight 

ID 

events  from  the  Fiji  Islands  were  normalized  using  (A-6)  with 
S  =  1  as  shown  in  Figure  A-l  The  average  gain  of  each  station 
was  computed  using  (A-22),  and  the  deviations  from  these  averages 
were  plotted  as  shown  in  Figure  A-2 .  This  data  was  then  broken 
into  class  intervals  and  plotted  on  a  logarithmic  probability 
graph,  Figure  A- 3. 


The  fit  to  a  log-normal  distribution  function  is  sufficiently 
good  to  justify  the  assumptions  that  the  additive  noise  is  neg¬ 
ligible  for  our  purposes  provided  we  use  large  events.  Also, 
since  the  error  term  in  the  linearized  model  is  distributed 
normally,  classical  least  squares  methods  with  the  usual  statis¬ 
tical  qualifications  on  all  estimates  can  be  applied  directly. 


Another  way  of  presenting  this  data  is  shown  in  Figure  A-4 . 

There  the  observed  amplitudes  for  these  events  are  plotted  against 

the  model  values,  namely  observed  log  L.  versus  log  S,  +  log  g  , 

lj  l  l 

where  S  and  g^  are  estimated  using  Equations  (A-20)  and  (A-21). 


A-7 


ssamndwv  Q3znvwaoN 


Figure  A-l  .  Eight  Fiji  Island  Events  -  Normalized  and  Showing  the 
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Figure  A-2 .  Eight  Fiji  Island  Events  Showing  Deviations  about  the 
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APPENDIX  B 


SEISMIC  DATA  LABORATORY 
ALEXANDRIA,  VIRGINIA 


DIGITAL  COMPUTING  SECTION 


A .  IDENTIFICATION 

Title ;  ANOMALY 

COOP  Identification;  Z108 

Category:  General 

Programmer :  Nicholas  Fletcher 

Date ;  16  June  1966 

B .  PURPOSE 

To  compute  travel-time  and  amplitude  anomalies 

C .  USAGE 

1.  Operational  Procedure:  This  is  a  FORTRAN  63  program. 
All  input  is  from  punched  cards  and  the  output  is 
from  the  printer,  and  an  off-line  card  image  tape. 

2.  Parameters  Input  parameters  consist  of  the  latitude 

and  longitude  of  the  location  of  the  subarray  centers. 

3.  Space  Required  Not  Applicable 

4.  Temporary  Storage  Requirements:  None 

5 .  Alarms ;  None 

6.  Error  Codes:  None 

7 .  Error  Stops:  None 

8.  Input  and  Output  Tape  Mounting:  There  is  no  input 
tape.  The  output  tape  is  on  logical  Unit  16. 

9.  Input  and  Output  Formats: 


B-l 


The  input  card  formats  are  as  follows: 


Subarray  locations  deck,  must  be  the  first  deck 
of  cards  read  by  the  program.  The  first  card  of 
this  deck  must  contain  the  number  of  subarray 
location  cards.  This  number  is  punched  in 
Columns  1  through  5  right  justified.  The  first 
card  is  followed  by  the  subarray  location  cards, 
with  one  location  per  card.  The  program  is 
limited  to  a  maximum  of  51  subarrays.  1  (One) 
punched  in  Column  10,  prints  out  J.  B.  Table  and 
station  locations. 

Subarray  location  card  format  is  as  follows : 

FORMAT  (A8,  2  (F5  0,  F3.0,  F5.1,  A1 ) 

COLS  . 

1-8  Subarray  name,  left  justified 

9-13  Subarray  latitude  degrees 

14-16  Subarray  latitude  minutes 

17-21  Subarray  latitude  seconds  (with  decimal  point 

punched ) 

22  N  or  S  for  north  or  south  latitude 

23-27  Subarray  longitude  degrees 

28-30  Subarray  longitude  minutes 

31-35  Subarray  longitude  seconds  (with  decimal  point 

punched ) 

36  E  or  W  for  east  or  west  longitude 

The  location  cards  are  followed  by  the  first 
event  card  - 

Event  Card  Format  is  as  follows : 

Columns  1  through  14  contain  the  origin  time  of 
the  event 


COLS  . 

1-2 


3-4 


Month 

Day 


B-2 


COLS  0 


5-6 

Year 

7-8 

Hour 

9-10 

Minute 

11-14 

Second  to  the  nearest  tenth  of  a  second 

15-18 

Latitude  of  the  event  in  degrees  to  the  nearest 
tenth 

19 

N  or  S  for  north  or  south 

2  0-24 

Longitude  of  the  event  in  degrees  to  the  nearest 
tenth 

25 

E  or  W  for  east  or  west  longitude 

26-28 

Depth  of  event  to  the  nearest  kilometer 

29-31 

Magnitude  to  the  nearest  tenth 

32-36 

This  is  an  output  field.  The  program  will  com¬ 
pute  the  distance  (to  the  nearest  kilometer) 
between  the  event  and  the  standard  subarray 
(AO  for  LASA) 

37-39 

This  is  an  output  field .  The  program  will  com¬ 
pute  the  station  to  epicenter  azimuth  (to  the 
nearest  degree  relative  to  AO) „ 

40-78 

Not  used  by  the  program,  but  reproduced  on 
output  0 

79-80 

Must  be  blank,  to  identify  this  card  as  an  event 
card  and  not  a  phase  card 0 

Each  event  card  should  be  followed  by  phase 
cards p  A  phase  card  contains  the  observed  phase 
data  for  a  given  subarray.  Each  phase  card 
must  have  an  index  number  punched  in  Columns 
of  the  subarray  location  cards  (See  Sec-C9) . 
There  does  not  have  to  be  a  phase  card  for 
every  subarray  location  card  nor  a  specific 
numerical  order* 

Phase  Card  Format : 

FORMAT  (Columns  1  through  9  contain  the  phase  arrival 
time  at  this  subarray) 
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COLS  . 
1-2 
3-4 
5-9 
10 

11-14 

15-16 

17-20 

21 

30 

31-35 

36-40 

41-42 

43 

44-45 

46-50 

51-55 

56-60 

61-78 

79-80 


Hour 

Minute 

Second  to  the  nearest  hundredth 
Not  Used 

Amplitude  (zero  to  peak)  to  the  nearest 
millimicron 

Not  Used 

Amplitude  (peak  to  peak)  to  the  nearest  millimicron 
Not  Used 
Signal  Quality 

Time  anomaly  to  the  nearest  hundredth  of  a 
second  —  this  is  an  output  item  computed  by 
the  program. 

Amplitude  anomaly  zero  to  peak  relative  to 
geometric  mean  computed  by  the  program  (output 
field) 

Subarray  name  code  (output  item) 

Not  Used 

Reference  subarray  (for  each  event)  (only  in 
last  phase  card) 

Peak  to  peak  amplitude  anomaly  relative  to 
geometric  mean  (output  field) 

Zero  to  peak  amplitude  anomaly  normalized  to 
average  amplitude  (0-P)  (output  field) 

Peak  to  peak  amplitude  anomaly  normalized  to 
average  amplitude  (P-P)  (output  field) 

Reserved  but  not  used  by  program 

Index  number  of  subarray  recording  this  phase 
must  correspond  to  subarray  location  card 

Additional  events  with  their  associated  phase 
cards  follow  in  like  fashion  for  as  many  events 
as  desired..  The  subarray  location  cards  are 
not  repeated  in  the  deck  setup. 
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The  output  tape  is  a  card  image  tape,  with 
the  same  information  and  format  as  the  input 
tape  with  the  addition  of  the  computed  values 
for  distance  and  azimuth  for  event  cards,  and 
travel-time  and  amplitude  anomalies  for  the 
phase  cards. 

1 0 .  Selective  Jump  and  Stop  Settings :  None 

11.  Timing  :  Approximately  3.5  seconds  per  event 

12.  Accuracy:  No  better  than  the  input  data 

13.  Cautions  to  User.  Every  phase  card  must  have  an  index 

number.  If  Columns  79-80  are  blank,  the  program 
assumes  it  is  an  event  card 


14  Equipment  Configuration.  Standard  F-63  COOP  monitor 
system  with  the  card  reader  on  Unit  7  and  the 
printer  on  Unit  6  and  the  output  tape  on  Unit  16. 

15.  References :  Letter  from  MIT  dated  April  8,  1966 


D .  METHOD 

Definitions 

Observed  Travel-Time  =  Observed  Phase  Arrival  Time  Minus 

Event  Origin  Time 

Predicted  Travel-Time  =  Travel-Time  for  the  J-B  Table. 

Found  by  table  look-up  on  distance  and  depth  with 
four-point  interpolation  for  distance  and  linear 
interpolation  for  depth. 

0  =  Observed  travel-time  at  subarray  i 

=  Predicted  travel-time  to  subarray  l 

0^_  —  Observed  travel-time  at  standard  subarray 

—  Predicted  travel-time  to  standard  subarray 

Transit  Time  Residual  (Travel-Time  Anomaly)  =0  -0  +H  -H. 

i  r  r  i 

(Relative)  Amplitude  Anomaly  =  Amplitude  (Zero  to  Peak)  at 
subarray  i  divided  by  geometric  mean  amplitude  (zero  to  peak) 
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(Average)  Amplitude  Anomaly  =  Amplitude  (Zero  to  Peak)  at  subarray  i 
divided  by  arithmetic  average  amplitude  (zero  to  peak)  for  all  funct¬ 
ioning  subarrays. 

Relative  and  average  amplitude  anomalies  are  also  computed  in  like 
manner  for  the  peak  to  peak  measurements. 


APPENDIX  C 


SEISMIC  DATA  LABORATORY 
ALEXANDRIA,  VIRGINIA 

DIGITAL  COMPUTING  SECTION 


A.  IDENTIFICATION 

Title ;  Display  and  Analysis  of  Anomalies 
COOP  Identification;  Z95  DISPLAY 
Category  General 
Programmer ;  J.  W.  Monroe 
Date;  19  July  1966 


B.  PURPOSE 

To  take  the  output  from  the  anomaly  program  and  com¬ 
pute  the  average,  standard  deviation,  and  number  of  occur¬ 
rences  of  a  selected  anomaly  for  selected  subarrays  or 
sensors,  distances,  and  azimuths. 

C .  USAGE 

1.  Operational  Procedure:  This  FORTRAN-63  program  is 
with  a  binary  or  symbolic  deck  and  a  set  of 
control  cards  for  each  use. 

2 .  Parameters  ^  None 

3.  Space  Required:  21,222  buffers 

4 .  Temporary  Storage  Requirement;  None 

5 „  Print-Outs :  The  program  first  prints  out  the  anomaly 
used,  followed  by  the  distance  and  azimuth  criteria 
and  the  subarrays  or  sensors  used.  If  any  seismo¬ 
grams  are  to  be  deleted  from  processing,  they  are 
listed  next.  Following  this  will  be  a  page  for 
each  subarray  or  sensor  with  the  average,  standard 
deviation  and  number  of  occurrences  of  the  selected 
anomaly  for  each  distance  and  azimuth  chosen.  If 
plotting  is  requested,  the  last  page  will  contain 
scale  factors  and  other  information  pertaining  to 
the  plots. 
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6 . 

Error 

Returns:  None 

7  . 

Error 

Stops:  None 

8. 

Input 

and  Output  Tape  Mountings: 

Logical  Unit  #4  -  Input  BCD  tape  from  ANOMALY 
Logical  Unit  #3  -  Output  plot  tape  if  desired 
9.  Input  and  Output  Formats* 

INPUT:  (only  one  case  per  run) 


Card 

#  Cols 

Format 

Descript  ion 

1 

1 

11 

Anomaly  to  be  used  (see  last  page) 

5 

11 

Plot  switch  1-plot,  blank-no  plot 

2 

1-10 

110 

Number  of  Dels  (Max  19) 

11-20 

110 

Number  of  Azimuths  (Max  13) 

21-30 

110 

Number  of  Subarrays  or  Sensors 

(Max  21) 

next  N 

1-10 

F10.2 

Distance  (km) 

11-20 

F10.2 

+_  Distance  (km) 

Where  N 

=  NO. 

Punched 

in  Cols  1-10  on  card  #2  (Max  19) 

next  M 

1-10 

F 10  „  2 

Azimuth  (Deg) 

11-20 

F10.2 

+_  Azimuth  (Deg) 

Where  M 

=  NO. 

punched 

in  Cols  11-20  on  card  #2  (Max  13) 

next  L 

1-10 

F10.0 

Subarray  or  Sensor  Number 

Where  L 

=  NO. 

punched 

in  Cols  21-30  on  card  #2  (Max  21) 

next  D 

1-4 

.14 

Seismogram  Nos.  to  be  deleted 

(Max  100) 

OUTPUT :  The  printed  output  will  be  that  described  in  PRINT-OUTS . 

The  plot  tape  on  logical  unit  #3  will  contain  a  plot  for 
each  subarray  or  sensor  of  the  average  for  that  subarray 
minus  the  average  of  averages  for  all  subarrays .  The  plotted 
values  will  be  computed  for  each  distance  at  each  azimuth 
(in  that  order)  until  N  points  have  been  plotted  where  N 
equals  the  number  of  dels  times  the  number  of  azimuths.  The 
plots  will  always  represent  a  new  del  every  tenth  of  an  inch. 
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10. 


D. 


Selective  Jump  and  Stop  Settings ;  None 

11.  Timings  2  milli-sec  per  case  where  each  case 

has  a  separate  distance,  azimuth,  subarray, 
and  anomaly. 

12.  Accuracy^  The  same  as  input  data 

13.  Cautions  to  Users;  Both  logical  units  3  and  4  must 

be  mentioned  in  the  COOP  card.  If  no  plot  is 
desired,  use  the  BY3  convention.  The  plot 
routine  is  not  internally  controlled  so  the 
user  must  request  the  same  number  of  files  to 
be  plotted  as  there  are  subarrays  or  sensors 
with  a  maximum  of  ten  across  the  plot.  Observe 
that  there  are  maximum  limits  on  the  number 
of  distance  and  azimuth  intervals  and  number 
of  subarrays  and  seismograms  to  be  deleted. 

14.  Equipment  Configuration;  Standard  COOP  for  FORTRAN-63 

15.  References  s 

Program  write-up  for  ANOMALY  Fletcher  16  June  1966 


METHOD 


Average  =  X  ,  =  - 

i , 3 ,k  N 


1*3 


Standard  Deviation  =  a 


Yx  k 

L  l , 3 ,k 

L  I  xi.i.lc  -CZxi.i.k)2/tli,i  ] 


i  j  .k 


N.  .  -1 


Number  of  Occurrences  =  N 

i ,  3  ,k 

where,  i  =  azimuth  +  deg  s 
j  =  distance  f  km’s 
k  =  subarray  or  sensor 

N  =  number  events  pertinent  to  the  specified  coordinates 
and  specified  subarray 

X  ,  =  anomaly  used 

l ,  3  ,k 


C-3 


SAMPLE  INPUT  DATA  FOR  DISPLAY 
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TO : 


Users  of  Program  DISPLAY 
FROM:  Jo  Wo  Monroe 

DATE:  2  September  1966 

SUBJECT:  Additional  Features  of  Program  DISPLAY 

REFERENCE:  Z95  DISPLAY  Program  Writeup,  19  July  1966 


The  user  may  now  run  DISPLAY  using  more  than  one  input  tape  de¬ 
scribed  in  the  program  writeup  referenced  above =  The  first  input  tape 
must  be  on  Logical  Unit  4  but  the  second  third  etc.,  tapes  may  be 
assigned  to  any  logical  unit  from  1  thru  49  inclusive,  except  Logical 
Unit  3  (the  plot  tape).  When  using  more  than  one  input  tape.  Logical 
Units  3,4,  and  all  other  logical  units  must  be  mentioned  in  the  COOP 
card.  The  control  card  setup  for  the  case  of  more  than  one  input  tape 
is  exactly  the  same  as  that  shown  in  the  SAMPLE  INPUT  DATA  FOR  DISPLAY 
of  the  writeup,  except  that  following  the  last  card  of  seismograms  to 
be  deleted  or  following  the  last  subarray  number  if  there  are  no  de¬ 
letion  cards,  a  blank  card  must  be  inserted  followed  by  the  logical 
tape  assignment  for  the  second  input  tape.  If  more  than  two  input 
tapes  are  used,  additional  cards  are  used  to  specify  the  logical  tape 
unit  for  input  tapes  3,4,  etc.,  but  the  blank  card  is  not  used  again. 
Format  for  these  tape  assignments  is  15.  The  following  is  a  sample 
card  setup  starting  with  the  deletion  cards: 


etc  „ 


BLANK  CARD 


^  1350 
(  4210 
^  3490 


fourth  input  tape 


_  logical  tape  unit  for  the 

third  input  tape 

logical  tape  unit  for  the  second 
input  tape 

signifies  the  end  of  the 
seismograms  to  be  deleted 
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As  soon  as  an  end  of  file  is  encountered  on  the  input  tape  being 
processed,  a  tape  assignment  card  is  read  in  and  the  program  pauses 
with  the  next  input  logical  unit  number  displayed  in  the  "A"  register. 
When  the  Start  Key  is  hit,  the  program  begins  processing  the  new  input 
tape.  Two  input  tape  units  might  be  used,  such  as  4  and  5  where  the 
first  tape  is  on  Unit  4  and  the  second  tape  is  mounted  on  5  while  4  is 
being  processed.  After  Unit  5  begins  processing,  another  tape  could  be 


mounted  on 

Unit  4. 

DEBUGGING  AIDS  FOR  DISPLAY 

Jump  Key  2 

On  displays  the  event  number  being  processed  in  octal 
in  the  "A"  register.  Program  halts  until  Start  Key 
is  hit.  For  exact  stops,  leave  JK-2  on  and  step 
through  events  with  Start  Key.  Program  halts  on  header 
card  of  event  being  displayed  in  "A"  register. 

Jump  Key  1 

Up  causes  printout  of  accumulation  that  far  into  the 
program.  Turn  JK-1  off  as  soon  as  printout  begins 

Jump  Key  3 

and  JK-3  off  as  soon  as  printout  has  finished  Use 

this  feature  sparingly  since  some  loss  of  accuracy 
occurs  when  sigma  is  set  equal  to  zero  rather  than  a 
negative  number.  All  other  times,  the  accumulation 
is  reset  to  its  original  value  without  loss  of  accu¬ 
racy. 

Be  sure  both  COOP  card  and  EXECUTE  card  have  sufficient  time  limits  for 
the  case  being  run.  For  a  case  of  324  events  with  21  stations,  17 
distances,  and  12  azimuths,  allow  at  least  60  minutes. 

CHANGE  IN  FORMAT 

The  format  of  the  control  cards  containing  the  numbers  of  seis¬ 
mograms  to  be  deleted  has  been  changed  from  14  to  15. 

I 
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APPENDIX  D 


The  product  error  model  described  in  Appendix  A  can  be 
used  to  detect  amplitude  data  which  is  inconsistent  with  the 
population  from  which  it  is  drawn,,  Several  approaches  to 
automating  this  error  detection  were  suggested  and  the  follow¬ 
ing  is  a  description  of  the  method  which  was  finally  selected 
for  production  usec 

The  collection  of  recorded  data  was  grouped  into  a  set 
of  cells .  The  grouping  geographically  was  by  distance  and 
azimuth  and  a  cell  is  defined  by  a  particular  distance  and  azi¬ 
muth  range  and  a  station  designation0  The  assumptions  made 
were  that  the  geographic  regions  chosen  are  sufficiently  small 
so  that  the  true  anomalies  are  constant  within  each  region 
and  that  the  variability  is  the  same  for  all  cells 

Using  the  results  of  Appendix  A,  namely  that 

2  ? 
log  e  ^  N  (0,  a  ),  one  may  make  estimates,  denoted  as  , 

2 

of  a  for  each  cello  These  can  be  pooled  to  obtain  a  very 

2 

good  estimate  of  a  ,  and  with  this  estimate  and  known  dis¬ 
tribution,  outlying  estimates  and  thus  probably  blunders  can 
be  detected, 

2 

The  estimates  of  a  ,  computed  for  each  cell  (given 
station  3  and  region)  are 


s2  - 
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where  I  is  the  number  of  events  recorded  in  the  particular  cell 
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X  =  log 
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2  2 
and  S  is  distributed  as  \ 
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By  pooling  these  variance  estimates,  a  very  good 
overall  estimate  is  obtained  due  to  the  large  number  of  data. 
The  estimates  in  each  geographical  cell  were  quantized, 
indexed  and  entered  on  a  histogram,  (Figure  6  of  the  main 
report)  as  a  function  of  the  number  of  events  in  that  cell. 

A  pooled  estimate  of  the  population  variance  is  then 
given  by 


s2 
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where  f  .  is  the  number  of  estimates  in  the  class  interval 
<  D 

indexed  by  j  and  with  I-events  in  the  cells.  The  value  of 
2 

S  is  equal  to  the  center  value  of  the  variance  in  the 


class  interval,  j. 

2 

Using  this  estimate  for  a  we  may  compute  100(l-a)% 

2  2  2 
confidence  intervals  for  S„  and  S  .  Thus,  we  set  a  =  S 

N  N 

and  define  confidence  intervals. 
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2 

n  \ 


—J 
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X  ( a/2 , f ) 


2 

X  (l-a/2,f) 
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or 
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x  (a/2, f) 

rr  <C  Q  ^ 

2 

x  (l-a/2, f 

f 

a  SN 

where  f  =  number  of  degrees  of  freedom,  (I  -  1), 
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( a/2 , f ) 


and 


2 

x  (l-a/2 , f ) 


are  found  in  tables. 
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